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The invaded cluster approach is extended to 2D Potts model with annealed vacancies by using 
the random-cluster representation. Geometrical arguments are used to propose the algorithm which 
converges to the tricritical point in the two-dimensional parameter space spanned by temperature 
and the chemical potential of vacancies. The tricritical point is identified as a simultaneous onset of 
the percolation of a Fortuin-Kasteleyn cluster and of a percolation of "geometrical disorder cluster" . 
The location of the tricritical point and the concentration of vacancies for q = 1, 2, 3 are found 
to be in good agreement with the best known results. Scaling properties of the percolating scaling 
cluster and related critical exponents are also presented. 



I. INTRODUCTION 

Detailed investigations of fractal properties related to 
criticality were done some time ago for the geometrical 
phase transitions such as percolation The interest 
for similar properties in the context of thermal phase 
transitions @, 0, less exploited until know, has been re- 
newed recently^, 0,0,0]. Monte Carlo (MC) studies of 
phase transitions in the last decade have given rise to sev- 
eral cluster algorithms, based on the Fortuin-Kasteleyn 
(FK) representation of the partition function ||. In ad- 
dition to their principal task to reduce the critical slow- 
ing down present in the local update algorithms, they 
have an advantage to offer a better insight into geomet- 
rical aspects of phase transitions and represent a natural 
tool for numerical studies of these phenomena. An al- 
gorithm in which this geometrical approach is used in 
a particular way is the invaded cluster (IC) algorithm, 
defined by Machta et al. @ . While in the standard clus- 
ter approaches, such as the Swendsen-Wang (SW) 
or Wolff algorithm |Tl| , the clusters are built for a given 
temperature, the IC algorithm starts from some geomet- 
rical property of the criticality that can be generated by 
a random process and obtains the critical temperature 
as an output. It was applied to Ising and Potts models 
and, later, in a series of other studies, e. g. on the fully 
frustrated Ising model [l2j , or the XY model [HI . It ap- 
pears equally efficient in both the second- and first-order 
phase transitions. 

The tricritical point present in systems which exhibit 
the changeover from the first- to second-order phase tran- 
sitions is difficult to access in numerical and finite-size 
scaling studies due to crossover effects. Even the loca- 
tion of the tricritical point appears to be a difficult task 
in many cases, from models with long-range interactions 
0, E Gil t° models with the quenched dilution [l7j]. 
Standard MC approaches identify this point as the onset 
of a first-order transition, recognized by two maxima in 
the free energy distribution which scale as surface [H, Il9| 



and can be analyzed directly, or from the Binder's fourth 
cumulant [20| . Some alternative numerical approaches 
were also proposed - e.g. the microcanonical MC study 
of the tricritical point [21j in Blume-Capel model. In the 
present work we examine a different approach and extend 
the IC algorithm by introducing an additional geometri- 
cal condition. We present it on the example of a 2D Potts 
model with the annealed dilution. 

The paper is organized as follows: In section|TT]we dis- 
cuss the model and its graphical expansion in the diluted 
case. In section lllll we present the extension of the IC ap- 
proach to the tricritical point and explain the algorithm. 
In section [Pvl we discuss our results for the location of the 
tricritical point, scaling properties of percolating cluster, 
and related critical exponents. Section V contains the 
conclusion. 



II. MODEL 

We consider the q state Potts model [22j with annealed 
vacancies on a square lattice, described by the Hamilto- 
nian 
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where s.j denotes the g-state Potts variable at the site i 
and the variable U takes the values or 1 when the site i 
is empty or occupied, respectively. The sum is taken over 
nearest neighbors, J > is a ferromagnetic coupling, and 
G is the chemical potential of vacancies. The pure model, 
where a single thermodynamic parameter - temperature 
- is governing a transition, is recovered in the limit G — » 
—oo. 

In the pure model the transition is of the second order 
for q < 4 and of the first order for q > 4 [23[ . The pres- 
ence of annealed vacancies induces the first-order phase 
transition, so that for q < 4 the transition line in the 
(T, G) parameter space contains both regimes of the first- 
and second-order phase transition, separated by a tricrit- 
ical point. 

For the pure 2D Potts model, the exact analytical ex- 
pressions exist both for the critical temperature and crit- 
ical exponents [2jJ. For a diluted case, exact analytical 
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expressions for tricritical exponents are also available, 
from the conformal theory [2 51 ] and Coulomb gas map- 
ping [26[, while for the location of the tricritical point 
there are only good approximate results [2^ |. 

As shown by Fortuin and Kasteleyn Q , the pure Potts 
model is equivalent to the random-cluster model, which 
may be understood as a generalized percolation model. 
The random-cluster partition function is given by 

z = J2p b{l) (i-p) E ~ bh) <i ch K (2) 

p=\-e~ 0J . (3) 

The summation runs over the set of all the graphs on the 
lattice r. Each graph represents one choice of placing 
the bonds on the lattice. The quantity p (where /3 is the 
Boltzmann factor) may be interpreted as the probability 
of presence of a bond on an edge, 6(7) is the number of 
bonds in the graph 7, E is the total number of edges 
(maximum number of bonds), which for the square lat- 
tice of the linear size L is equal to 2L 2 . The entropy 
factor for each cluster is given by q, and 0(7) is the num- 
ber of connected components of a graph. FK clusters 
have a physical meaning - the probability that two spins 
separated by a distance r are in the same FK cluster is 
proportional to the correlation function. 

Graph expansion for diluted Potts model has already 
been studied elsewhere [28|, [29[ with a different choice of 
a Hamiltonian. It is easy to perform a graph expansion 
for the Hamiltonian fT]) by proceeding along the same 
lines as in the pure model. By using the equality 

JJi^Uti = [1 + {e PJ _ 1)Ssi ^ t . tj]t (4) 

the partition function can be expanded over all the pos- 
sible graphs on the lattice 

^ = ^ J EEEII e " /3G(tl " 1) II ^-i)s susMj , 

ier{si}{ti} i <i,j> 

(5) 

where {si} and {ti} denote the summation over all the 
configurations of Si and ti respectively. The structure 
of graphs remains the same as in the pure model. In 
each graph 7 one may separate the summation over the 
configurations including sites belonging to the clusters of 
sizes > 1 from the summation over the "isolated sites" 
consisting of single spins and vacancies. The summation 
over isolated sites may be seen as a lattice gas of single 
spins and vacancies in the field and gives 

E M ^""g" """" = (e? G + q) n " , (6) 

n v — ^ v ^ 

where 71,3(7) is the number of isolated sites on the graph, 
and n v is the number of vacancies. The partition function 
may be written in a more condensed form 

Z = £y(7)(l _p)S-6(7) g £(7) (e /3G + g) n ia ( 7 ) j (7) 

7er 



to be compared with the pure case Eq. ((2|). In the first 
three factors we recover the same expression as for the 
pure model except that 0(7) denotes only the number 
of connected clusters and excludes single spins. In the 
limit G — > — 00 Eq. ([7]) reduces to Eq. J2]) for the pure 
case. The only change occurs in the contribution to the 
entropy factor stemming from the isolated sites in each 
graph. While in the pure case it was equal to q nis , in 
the presence of dilution it corresponds to the one of a 
lattice gas, given by Eq. ©. The chemical potential 
G, responsible for dilution can be expressed through the 
parameter 



which, according to Eq. ([6]), has the meaning of the a 
priori probability to find a single spin on an isolated site 
on a graph. 



III. METHOD 

In Swendsen-Wang and Wolff algorithms one uses the 
knowledge of temperature to construct FK clusters by 
adding bonds with probability p. Clusters can then be 
independently flipped to another Potts state. The oppo- 
site happens in the IC algorithm of Machta et al. FK 
clusters are grown by a procedure inspired by the inva- 
sion percolation (30j . Bonds are placed at random be- 
tween neighboring spins in equal states, with probability 
1, until a geometrical condition, imposed on FK clusters 
(e.g. the bond percolation), is achieved. The tempera- 
ture is then deduced by equating the bond probability p 
defined in Eq. ([3]) with the ratio of the number of bonds 
to the number of satisfied neighbors. In particular, if the 
geometrical condition is the onset of the percolating FK 
cluster, the output temperature converges to the tran- 
sition temperature in the thermodynamic limit. In the 
language of connectivity the appearance of the percolat- 
ing FK cluster on the lattice corresponds to the state 
when the correlation length reaches the system size. 

The tricritical point is determined by two parameters, 
the tricritical temperature T t and field Gt, so the algo- 
rithm must include an additional condition besides the 
percolation of an FK cluster in order to converge to it. 
In the present 2D case we set the additional condition to 
be the percolation threshold of a "geometrical disorder 
cluster" defined as a cluster consisting of vacancies and 
single spins. That cluster is geometrical in the sense that 
two adjacent sites containing either a vacancy or a single 
spin are considered to belong to the same disorder clus- 
ter with probability 1 . This choice of the disorder cluster 
is justified because single spins behave in the same way 
as vacancies as far as correlations are concerned. It is 
supported by the picture obtained from the real-space 
renormalization group where the single spins take part 
in the renormalization of vacancies 13111. 
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One can also explain the second condition for the tri- 
critical point in terms of the persistence length £ [H, [33| . 
At the line of second-order transitions the correlation 
length £ diverges, but in the case of a tricritical point, 
there is another length scale in the system, £, which is 
related to the size of disorder clusters [H] . The length £ 
is finite on the line of second-order transitions but it di- 
verges at the tricritical point and remains infinite on the 
line of first-order transitions. Our algorithm generates 
finite system configurations in which £ oc £ oc L, since 
the percolation of a geometrical disorder cluster can be 
interpreted as £ reaching the size L. 



A. Algorithm 

An algorithm that is meant for locating the tricriti- 
cal point should find the percolation of the FK and ge- 
ometrical disorder clusters at the same time. In this re- 
spect we notice that in 2D the simultaneous percolation 
of two clusters is topologically obstructed unless it occurs 
in only one direction. Such anisotropic cases occur in a 
finite lattice and will be facilitated by a diverging persis- 
tence length, which occurs when approaching a first-order 
transition regime. 

Let us outline the basic steps to be explained in detail 
below. We start the MC iterations with random configu- 
ration of spins and randomly distributed vacancies with 
some initial concentration. Each MC iteration includes 
three steps: (i) the formation the FK clusters in the same 
way as the IC algorithm of Machta et al. H, (ii) iden- 
tification of disordered clusters including a check for the 
percolation, (iii) randomization of FK clusters and va- 
cancies (keeping the concentration unchanged). The re- 
moving or adding of vacancies, necessary to achieve the 
percolation threshold of the geometrical disorder cluster, 
is performed in more spaced intervals of MC iterations, 
and by imposing limitations on the number of vacancies 
created (destroyed) as described in Sect. IIII A 41 Namely, 
when the number of spins does not change, both the mean 
temperature and the field self-regulate to the transition- 
line. By changing the number of single spins sufficiently 
slowly it can be driven along the transition-line to the 
percolation threshold of the geometrical disorder cluster. 



1. Formation of FK clusters 

The simulation starts with some configuration of the 
Potts spins and vacancies on the lattice. Bonds are 
placed randomly, between the Potts spins in equal states. 
When percolation is achieved, or all the edges have been 
probed, the procedure is terminated giving a graph con- 
sisting of FK clusters. If percolation is achieved, this 
means that a quasi-critical configuration has been found 
for a given number of vacancies. 

There are several ways to characterize the onset of a 
percolating cluster during simulations. Two examples 



are: (a) winding of a cluster around a lattice with peri- 
odic boundary conditions; (b) the span of a cluster be- 
comes equal to the lattice size. The first characteriza- 
tion, termed a topological percolation, will be used in 
this work - from now on, when percolation on a finite 
lattice is mentioned it is meant in the sense of (a). 

The percolation of an FK cluster is determined by a 
procedure which uses the connecting vectors, described 
inRcf. [|. 

2. Identification of geometrical disorder clusters 

Identification of geometrical disorder clusters is similar 
to the formation of FK clusters. The edges are probed 
and adjacent isolated sites are connected by imaginary 
bonds in order to identify all the geometrical disorder 
clusters. The disorder bonds are said to be imaginary 
in the sense that they do not have a physical meaning. 
The formation of geometrical disorder clusters terminates 
when all the edges have been probed. Percolation is de- 
tected in the same way as for the FK clusters. 

3. Randomization 

Regardless whether the quasi-tricritical configuration 
was found or not, after the two cluster formation steps we 
randomize the configuration and iterate the procedure. 
The FK clusters in the configuration are randomized by 
flipping each of them into a new randomly chosen Potts 
state. The single spins and vacancies are randomized 
by exchanging their positions with probability 1. This 
procedure conserves the number of spins in the system. 

4- Adding (removing) vacancies and thermalization 

As far as temperature is concerned, the IC algorithm is 
self-regulating, as is easily understood in the pure model. 
If it starts with a configuration corresponding to T < Tc, 
the FK clusters are expected to percolate after a rela- 
tively small amount of bonds has been placed, because 
there is a small number of obstacles for the formation 
of a percolating cluster. They are caused by the bound- 
aries between the geometrical Potts clusters. The config- 
uration thus generated will correspond to a temperature 
higher than the starting one. When applied to a config- 
uration corresponding to T > Tc, with small clusters, 
in different Potts states, the IC algorithm meets a lot of 
obstacles. The number of bonds that need to be placed 
is quite large, corresponding to a temperature which is 
lower than the starting one. 

Despite the additional obstacles caused by vacancies, 
the temperature self-regulation remains in the present 
algorithm. However, we were not able to find a simi- 
larly elegant way to self-regulate the formation of the 
disorder cluster. Namely, the fact that the concentration 
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of vacancies at the tricritical point is far below the site 
percolation threshold, indicates that the percolation of 
geometrical disorder cluster is not random, but strongly 
correlated. Consequently, the thermalization procedure 
is important. In order to keep the number of vacancies 
just at the minimum necessary for the percolation of the 
geometrical disorder cluster, the system has to be allowed 
to relax each time that number of vacancies was changed. 

To this purpose we introduce the intervals of r MC 
steps during which the concentration of vacancies may 
vary by limited, very small amount. For each interval r 
we record the number of times b that the geometrical dis- 
order cluster has percolated together with an FK cluster. 
Further, we set a condition on the fraction b of simulta- 
neous percolation events during the r MC steps, to be 
less than certain small value bo. To obtain the threshold 
of disorder cluster percolation, the fraction bo needs to 
tend to 0, but for numerical reasons we set it to a small 
value of the order of 10~ 3 to 10~ 2 . 

If b > bo after r MC steps, a small, limited number 
of vacancies is removed and replaced by spins during the 
next r steps. If 6 is after r MC steps the same limited 
number of single spins is switched into vacancies during 
this period. If b < bo, the number of spins remains un- 
changed for the next r MC steps. The number of spins 
that is added or removed has to be small enough not to 
disturb the effective correlation established between va- 
cancies in geometrical disorder clusters of the tricritical 
point. 



5. Geometrical parameters 

When quasi-tricritical condition is fulfilled, all the rel- 
evant data are recorded and used for statistics; if quasi- 
tricritical configuration has not been found, no record is 
done. 

The thermodynamic parameters T and G are calcu- 
lated from the geometrical quantities. As in the pure 
case, the temperature is calculated from the bond prob- 
ability defined in Eq. ([3]), which can be expressed by the 
ratio 



P 



(9) 



where rif, denotes the number of bonds on the lattice and 
n ss the number of neighboring pairs in the same state. 

The chemical potential G is expressed through the a 
priori probability z defined by Eq. ([5]). It can thus also 
be expressed as the average fraction of single spins (rid) in 
the "lattice gas" of isolated sites consisting of vacancies 
and single spins, 



rid 



rid 



(10) 



We calculate the averages < p > and < z > , where < > 
is taken over the ensemble generated by the algorithm. 



The distribution of variables generated by the algorithm 
is not canonical, but is sharply peaked around tricritical 
values since it has a feedback mechanism that adjusts 
both parameters to the tricritical point. This was already 
pointed out earlier 0, [34| in the case of IC algorithm for 
the pure Potts case. 



IV. RESULTS 

We have considered the dilute Potts model described 
by the Hamiltonian ([1]) in three different cases, q = 1, 2, 
and 3. 

The calculations included lattices of linear sizes 24 < 
L < 240 for the cases q = 1 and 2, and 24 < L < 120 for 
the case q = 3. The statistics used was 3 • 10 5 MC steps 
for smaller and up to 5 • 10 4 MC steps for larger lattices. 
The statistics was also reduced to 10 4 MC steps in the 
case q = 3. 

The intervals r described in section III. A. 3. are cho- 
sen to be from 5 • 10 2 to 10 3 MC steps and the allowed 
fraction bo of the order of 10~ 2 . This value of bo was 
found sufficiently small not to produce systematic shift 
towards the first-order transition which would exceed sta- 
tistical error bars. 

The number of vacancies allowed to be changed after r 
MC steps is chosen to be L/24 for q = 1 and 2 and L/12 
for q = 3. The choice is arbitrary but allows the runs 
to be made on smaller lattices such as defined above and 
give reasonable accuracy. Also notice that the ratio of 
the number of single spins to be changed to the average 
number of single spins in the system tends to zero as l/L, 
so the non-equilibrium effects decrease as L increases. 

We present first the finite-size results and their extrap- 
olations for the concentration < t > and parameters p 
and z at the tricritical point. Further, we give the anal- 
ysis of the scaling properties related to the percolating 
clusters and evaluate several tricritical exponents. The 
error bars of Monte Carlo statistics in presented figures 
are less than symbol sizes. 

The results for tricritical exponents are compared to 
the exact values known from the conformal theory [25|, [35[ 
and Coulomb gas map [26], while the tricritical values 
of parameters are compared to very accurate results ob- 
tained by the transfer matrix technique [|| . 



A. The concentration 

The quantity for which we have obtained the best pre- 
cision is the tricritical concentration of spins, < t > t . 
The finite-size values < t >l arc displayed in Fig. (fTJ) as 
a function of the system size L. 

Simple fits to the power-law form < t >l=< t >t 
+b ■ L~ x were sufficient to obtain the extrapolations to 
L — > oo, which agree to the third digit with the best 
known results, as shown on Table I. 
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Figure 1: (Color online) Plot of < t >l vs. 1/L for different 
values of q. Plain lines denote the nonlinear fits to the power 
law form < t >l=< t > t +b ■ L~ x . 



Table I: Extrapolated values obtained for the tricritical concentra- 
tion 



q 


< t > 


< * > a 


l 


0.499 ± 0.001 


1/2 


2 


0.579 ± 0.002 


0.57979(1) 


3 


0.652 ±0.005 


0.65423(4) 



"Best known results from Ref. [y| 



B. Tricritical parameters 

The temperature and chemical potential are expressed 
here in terms of the geometrical parameters, the proba- 
bilities p and z of the random-cluster expansion i[T]). and 
are calculated using Eqs. (|9|) and (flUf . 

In Figs. ([2]) and ([3]) we display all the three sets of 
finite-size results for p and z, respectively. They are pre- 
sented in the form of difference of finite-size results and 
the expected tricritical values calculated using the pa- 
rameters given in Ref. [||. 

The accuracy of any extrapolations is lower by an or- 
der of magnitude compared to the one for the particle 
concentration. The reason for larger discrepancy may 
be attributed to the fact that these parameters exhibit 
much larger oscillations than concentration does during 
simulations and are more susceptible to the thcrmaliza- 
tion effects. This can be improved by a simple increase 
of statistics. 

As it may be observed in Figs. (J5J) and ([3]), the con- 
vergence of parameters pl and zl can be approximated 
rather well by a simple power-law form, which can be 



related to one of the critical exponents. 

The results for the convergence exponent x p , defined as 
the leading correction for the parameter p t (p(L) — p t oc 
L~ Xp ) are presented in Table II for different number of 
states q. They agree well with the values of the sublead- 
ing critical exponent y\. 

Notice that, in contrast to the ordinary critical point, 
where the convergence of the critical parameter at criti- 
cality is governed by the only relevant critical exponent 
yi , in the tricritical case there are two relevant exponents 
which both contribute to finite-size corrections of the crit- 
ical parameters. In the asymptotic regime the dominant 
contribution will come from the subleading one and not 
from y \ . 



Table II: Obtained values for the subleading tricritical exponent 
y% compared with exact values 









q 


Xp 


V2 


l 


1.02 ±0.05 


1 


2 


0.69 ±0.07 


4 

r, 


3 


0.54 ±0.03 


4 
7 



"Exact results from Ref. 



C. Scaling of FK clusters 

The fractal dimension of the percolating FK cluster is 
directly related to the anomalous dimension of the order 




« i i i i i i i i i i i i i 

3 3.5 4 4.5 5 5.5 6 
In L 

Figure 2: (Color online) Double logarithmic plot of | p(L) — 
pt | vs. L for different values of q. The best known values for 
Tt are given in Ref. [f| and values pt can be calculated to be 
0.8284271, 0.820168(1) and 0.807931(6) for q = 1, 2, and 3, 
respectively. 
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Figure 3: (Color online) Double logarithmic plot of | z(L)—zt 
vs. L for different values of q. The best known values for Gt 
are given in Ref. @] and values z t can be calculated to be 
0.02859548, 0.064220(5) and 0.113695(5) for q = 1, 2, and 3, 
respectively. 



parameter and equal to the magnetic critical exponent 
Uh- Its value at the tricritical point should be different 
from the one at the second-order transition line, while in 
the first-order transition regime it is simply equal to the 
embedding dimension. 

In Fig. |4]we present the log-log plots of the mass of the 
percolating FK clusters versus the lattice size. In Table 
III we present the obtained values for compared to 
the exact critical and tricritical values. However, in all 
the three cases considered here, the differences between 
critical and tricritical values arc rather small and of the 
order of the error bars of our calculations. We can only 
remark, that the results show a tendency to converge 
to higher values when only higher sizes are taken into 
account, which suggests that the considered percolation 
clusters are indeed tricritical. 

More conclusive information, whether the point in the 
parameter space produced by our algorithm corresponds 
indeed to the tricritical point and not to a point on the 
second-order transition line, can be obtained by examin- 
ing the scaling of the red bonds of the percolating FK 
cluster. 



D. Red bonds 



Figure 4: (Color online) Double logarithmic plot of perco- 
lating FK cluster size s max (L) vs. L for different values of 

q- 



Table III: Obtained values for the magnetic exponent compared 
with exact values for the critical and tricritical cases. 



q 


Vh 


Vh 




Vh 


l 


1.89 ±0.01 


1.89 ±0.02 


91 

48 


187 
96 


2 


1.88 ±0.02 


1.90 ±0.02 


15 
8 


77 

40 


3 


1.84 ±0.03 


1.88 ±0.02 


28 
15 


80 

42 



"Magnetic exponent obtained with all values of L 
''Magnetic exponent obtained with L > 72 for q - 



1 and 2 and 



L > 48 for q = 3 
c Exact results from Ref. 



3] 



they scale with a negative exponent, announcing the for- 
mation of compact clusters in the neighboring first-order 
transition regime. The red bonds corresponding to the 
ordinary critical point of the Hamiltonian {1} are known 
to have fractal dimensions of 3/4, 13/24, and 7/20 for 
q=l,2, and 3 respectively. 

Our finite-size results presented in Fig. O when fitted 
to a simple power-law form, give the values cited in Table 
IV. 

As it should be expected for relatively small sizes con- 
sidered in this study, the finite-size effects are still too 
strong to allow the calculation of negative exponents with 
sufficient precision, but the results clearly show negative 
exponents in all the tree cases. 



Red bonds are defined as the bonds which, when bro- 
ken, divide the whole cluster to which they belong into 
two parts [3(| . While at the second-order transition point 
the red bonds of the percolating FK cluster have some 
nontrivial fractal dimensionality, at the tricritical point 



V. CONCLUSION 

We have extended the invaded cluster algorithm to the 
calculation of the tricritical point. The algorithm is self- 
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Figure 5: (Color online) Number of red bonds on the perco- 
lating FK cluster vs. L 

Table IV: The red bond exponent at the tricritical point compared 
with exact values 



q 


XRB a 


Xrb 


X RB 


l 


-0.40 ± 0.06 


-0.50 ±0.08 


5 
8 


2 


-0.47 ±0.04 


-0.51 ±0.07 


19 

40 


3 


-0.11 ±0.07 


-0.20 ±0.08 


9 
28 



"Extrapolated form the complete set of values L 

'Extrapolated form values L > 72 for q = 1 and 2 and L > 48 for 



adjusting and locates the position of the tricritical point 
as the point of a simultaneous onset of the two percolat- 
ing clusters describing the percolation of order and perco- 
lation of disorder due to vacancies. While in the temper- 
ature variable the algorithm is completely self-adjusting, 
like for the simple criticality, in the parameter conjugate 
to the vacancy concentration the algorithm needs addi- 
tional fine tuning, in order to insure that the percolation 
threshold of the disorder cluster is achieved with mini- 
mum concentration of vacancies. 

The algorithm was illustrated on the example of a di- 
lute Potts model for three different values of q. It pro- 
duced the tricritical concentration with a good precision, 
and gave reasonably accurate results for the two tricriti- 
cal parameters. The analysis of the scaling properties of 
the obtained percolating clusters showed the characteris- 
tics proper to the tricritical point. The critical exponents 
were found in agreement with the tricritical exponents of 
the considered model. 

Let us mention, in the end, that the presented geomet- 
rical condition is not the only stopping rule within this 
invaded cluster procedure that might be constructed for 
locating the tricritical point. In future, it would be of 
interest to examine other possibilities (involving e.g. the 
red bonds), which could be applicable equally to higher 
dimensions, where the geometrical condition presented 
here, for topological reasons would not apply. 



q = 3 

c Exact results Ref. [4J 
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